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In the current density functional theory of linear and nonlinear time-dependent phenomena, the 
treatment of exchange and correlation beyond the level of the adiabatic local density approximation 
is shown to lead to the appearance of viscoelastic stresses in the electron fluid. Complex and 
frequency- dependent viscosity /elasticity coefficients are microscopically derived and expressed in 
terms of properties of the homogeneous electron gas. As a first consequence of this formalism, we 
provide an explicit formula for the linewidths of collective excitations in electronic systems. 
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Time-dependent density functional theory (TDFT) (!]] 
is frequently invoked as a tool for studying the dynamics 
of many-particle systems. This theory maps the difficult 
problem of interacting electrons in a time-dependent ex- 
ternal potential V(r, t) to the simpler one of noninter- 
acting electrons in an effective time-dependent potential 
v eff(r,t) = V(f,t) + v H (f,t) + v xc {r,t) (where v H is 
the Hartree potential, and v xc is the exchange-correlation 
(xc) potential) yielding the same density n(r,t). In or- 
der to obtain a practical computational scheme, the xc 
potential is usually approximated as a function of the 
instantaneous local density, 



ALDA 



(n{r,t)) = 



de xc {n) 
dn 



(1) 



n—n{r,t) 



where e xc (ri) is the xc energy density of the homogeneous 
electron gas of density n. This scheme is known as the 
"adiabatic local density approximation" (ALDA) M . Be- 
cause the conditions of validity of the ALDA (slowly vary- 
ing density and potential in time as well as in space) are 
seldom met in experiments, a few attempts have been 
made J3|,^| to improve upon the ALDA. The objective of 
these attempts was to obtain approximations for the xc 
potential which would still be local in space, but not in 
time. All these approximations were found to suffer from 
inconsistencies, such as the failure to satisfy the so called 
"harmonic potential theorem" (HPT) Q-^j], or other ba- 
sic symmetries. Only recently, it has become clear that 
the root of these difficulties lies in the fact that the xc 
potential in TDFT is an intrinsically nonlocal functional 
of the density, that is, a functional that does not admit 
a gradient expansion . 

Fortunately, a local gradient expansion is still possible 
if the theory is formulated in terms of the current den- 
sity. This was shown in ref. by Vignale and Kohn, 
who developed the time-dependent current density func- 
tional approach to the linear response theory, and gave an 



explicit expression for the linearized xc vector potential 
o lc (r, to) for a system of slowly varying density, subject 
to a spatially slowly varying external vector potential at 
a finite frequency u>. Their expression becomes exact in 
the limit k <C tu/vp,kp and q <C Lo/vF,kp, where fc _1 
and g _1 are the characteristic length scales for variation 
of the external potential and equilibrium density, respec- 
tively, and kp and vp are the local Fermi momentum 
and velocity. However, the final expression for a xc (r, u>) 
in ref. is rather formidable, and its physical meaning 
is far from transparent. Furthermore, it is restricted to 
the linear response regime. It is the purpose of this paper 
to overcome these limitations. 

In the following we derive a consistent local theory 
of the nonlinear dynamical response of a quantum elec- 
tronic system of "slowly varying density" , in the sense 
specified above. The effect of the xc potential beyond 
the ALDA will be shown to be analogous to the intro- 
duction of viscoelastic stresses in classical fluid dynamics 
and elasticity theory ||. The generalized viscosity coef- 
ficients (or, equivalently, the generalized bulk and shear 
moduli) are complex and frequency-dependent functions 
of the density, and can be calculated in terms of the lo- 
cal field factors of the uniform electron gas jl(J[fl). An 
important consequence of the xc viscosity is to provide a 
damping mechanism for long-lived collective excitations 
which cannot effectively decay into particle-hole pairs - 
the only form of damping allowed within the ALDA jl2) . 

We begin by recasting the linear response theory of 
|t],H in a form that is suitable for the nonlinear gener- 
alization. Let no(r) be the ground-state density of the 
system, and let XKSjijir,? 1 ,u>) be the current-current re- 
sponse function of a system of noninteracting electrons 
whose ground-state density is also no(r). The system is 
perturbed by a weak external vector potential d\ (r , t) = 
ai {f, bj)e~ luJt [jl3| , and we want to calculate the amplitude 
of the current-density response j\ (r , t) = ji (f, u})e~ iuJt , to 
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first order in a±. The answer is 



3 

+OHij(f',w) + a xc ij(r, uj)] dr, 



(2) 



where a#i is the first-order change in the Hartree po- 
tential (written in vector-potential form [fL3|), and a xc \ 
is the first-order xc vector potential, which contains the 
many-body effects. Note that a xc i in general has both 
longitudinal and transverse components, even when a\ 
is purely longitudinal. The local-density approximation 
for a xc \ was first derived in ref. @, eq. (19). We have 
found that that complicated formula can be written in a 
physically transparent form. To this end, we introduce 
the xc "electric field" E xc i(r,uj) = ^a xc i(r, uj). Then 
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7i (r) 4^ 



c9r 



(3) 



(e is the absolute value of the electron charge, and c is 
the speed of light). The first term is the linearization of 
the ALDA expression ([!]), and the dynamical correction 
is the divergence of the visco-elastic stress tensor 



dui duj 



-V • udi 



( XC V • uSi 



(4) 



Here u(r,uj) = ji(r,u))/rio(r) is the velocity field, and 
fj xc (uj, no(r)) and Cxc{u,no(r)) are complex viscosity co- 
efficients. They are related to the homogeneous electron 

gas functions f x l c L( UJ -> n ) an< ^ ■/•re t( w ' n ) (L f° r longitu- 
dinal, T for transverse) used in ref. as follows: 



( xc (u,n) =- — 

IUJ 



rh _ ^ rh _ ^ 

J xc.L xc,T dn 2 



(n) 



and 



■q xc {uj,n) = -—f 

IU1 



tc,T ■ 



(5) 



(6) 



The functions f£ c L r T \(uJ, n) are defined in terms of 
the dynamical local field factors C?l(t)(&, w) Q as 
fxc,L(T)( UJ ^ n ) = -lim fe ^ 47re 2 G L(T) (fc,w)/fc 2 . They 
have been recently calculated, within a mode-coupling 
approximation scheme, at zero temperature, by Conti 
et al. jlO],[ll]]. Because the local field factors are singu- 
lar for small k and u>, it is important to keep track of 
the order of the limits k — > and ui — » 0. The f xc 's 
are defined by taking the limit k — » first. Thus the 
limit for uj — » of f' x l cL (uj,n) differs from the famil- 
iar limfc^o 47re 2 Gi(fc, 0)/fc 2 = —d 2 e xc (n)/dn 2 . However, 
rigorous low frequency limits jl5| can be obtained 



from an analysis of the transport equation in the Landau 
theory of Fermi liquids [jl7| . We find 



lim Re 



and 



eh i \ 4 eh i \ d 2 e xc (n) 
/^(w,n) - -^fx C ,T(uin) j-±-L 



lim Re/* 

uj— >0 



2£ F F 2 /5 - F\J3 
5n l + Fx/3 ' 



(7) 



(8) 



where Fq,Fi... are the usual dimensionless Landau pa- 
rameters of the homogeneous electron gas |17[| . As for 
the imaginary parts, one finds Im/£ c L ^ T j ~ — c o,L(T) a -' 
for uj — ► 0, where the approximate values of the coeffi- 
cients Cq,l(t) are tabulated in ref. Q. 

From these results, one concludes that the real parts 
of the xc viscosity coefficients (which agree with the or- 
dinary notion of fluid viscosities) have finite values in 
the limit uj — ► 0. The imaginary parts of the viscos- 
ity coefficients are better understood in terms of bulk 
and shear moduli of an isotropic elastic medium, K~f{ n 
and ^ x ^ n . According to elasticity theory ||, we define 
Kll n {uj) = ujlmC xc and = wlm7fc c . The su- 

perscript dyn is a reminder that these are dynamical 
contributions to be added to the usual static ones, al- 
ready present in the ALDA. The static elastic constants 
are K^" 1 = n 2 d 2 e xc (n)/dn 2 and /4' c at = respectively. 
Equations (f?],||) show that for uj — > K^ n vanishes, while 
/i^™ has a finite value. A similar state of affairs holds 
for the noninteracting kinetic contributions to the bulk 



and shear moduli: K 



dyn 



and Mfein 



p(n), where 



p(n) is the noninteracting Fermi pressure (see below), 
and Kfcfo — ndp(n)/dn, = 0. The general conclu- 

sion is that dynamical (post- ALDA) effects do not mod- 
ify the bulk modulus, but they cause the appearance of 
a nonvanishing shear modulus and viscosity. 

Equations (Q) and ^ clearly display the basic symme- 
tries which were used in the derivation of First of 
all, the fact that the force exerted by the xc potential per 
unit volume — eno(r)E xc i(r, uj) — 711(7% uj)\7v xc o(r) can be 
written as the divergence of a symmetric rank 2 local 
tensor guarantees that the net force and the net torque 
acting on a volume element of the fluid have no contribu- 
tion from the volume element itself (Newton's third law) . 
Therefore, the HPT, the "zero force" and "zero-torque" 
theorems of [^J|] are manifestly satisfied. (Note that the 
force exerted by the xc "magnetic field" B xc = V X Q>xc 1^ 
legitimately disregarded in this argument, being of higher 
order in the gradient expansion [see also below]. Besides, 
it is rigorously absent in the linear response theory, if 
there is no static magnetic field). 

Thus far, we have used the condition of slow density 
variation (k,q <C uj/vf) only in approximating the xc 
vector potential. If this condition is met in the physical 
system under study, we can also use it to approximate 
the Kohn-Sham (KS) response function, XKS,ij{r, f, uj). 
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Then eq. (0) for the current reduces (after considerable 
algebraic manipulation) to a linearized Navier-Stokes 
equation |j] with complex and frequency-dependent vis- 
cosity coefficients: 

r c 

- imu)jn(r, lu) = n (r) -zw-a^r, uj) 
L c 



n 

■E 

3 



d<Jij{r,u>) 



(9) 



Here pi(n,) = p(n) — p(no) is the first order change in 
the pressure of the noninteracting electron fluid (p(n) = 
(3Tr 2 ) 2 ^h 2 n 5 ^/5m), and v H1 , V^ DA are the first-order 
changes in the Hartree potential and ALDA xc potential. 
The full stress tensor cr^ is defined as in eq. (|J) , with the 
viscosity coefficient fj xc replaced by f\ xc — p(no)/iu), while 
C, xc remains unchanged. 

If, on the other hand, the conditions k,q <C uj/vf are 
not well satisfied by our system, then it is better to revert 
to the original KS formulation (||) , which treats the non- 
interacting response exactly. The use of the local density 
approximation (||) for the xc potential becomes then an 
uncontrolled approximation, but it may still work well in 
practice. In particular, we note that eq. (|^) allows for the 
phenomenon of Landau damping (damping of collective 
modes by single electron- hole pairs), while eq. (^) does 
not. 

Let us now discuss the generalization of the formalism 
to the nonlinear response regime. In this case, one must 
solve the full time-dependent KS equations for the KS 
orbitals ip a (r , t) 



ifi 



0_ 

dt 2m 



1 / -» & e \ 2 

— \^-ihV + -<z(r, t) + -a xc (r, t) 



-v (r) - v H (r,t) 



l/> a (?,t)=0, (10) 



starting, for example, with the static KS orbitals cor- 
responding to the external potential vo(r) at the initial 
time. The density and the current density are computed 
from the KS orbitals according to the usual rules || . The 
form of the nonlinear xc vector potential is dictated by 
the following requirements: (i) the xc force density 



d _ - 
— + u • V | << , , . , 
at 



- 22 MfcVjQxcy 



(11) 



must be the divergence of a local symmetric rank two 
stress tensor (Newton's third law). 'Locality' here means 
that <T XCt ij(f, t) is a function of n(R, t'), j(R, t') and their 
spatial derivatives, where t' < t, and R(t'\r, t) is the posi- 
tion at time t' of the fluid element which evolves into r at 
time t Jl8| ]. (ii) Under transformation to an accelerated 
frame of reference [|j with origin at x(t), the stress tensor 



(r,t) becomes a' xclJ (f,t) {<r + x(t),t). (iii) 

Equation (|l^) must reduce to the linear response theory 
in the limit of small external perturbations, and to the 
nonlinear Navier-Stokes equation in the limit of slowly 
varying (in time) perturbation Jl9[ . 

To within the accuracy of our approximation, i.e., to 
second order in the spatial derivatives, the above require- 
ments uniquely determine the form of a xc : 



e da XCt j(r, t) _ alda 

c dt - VlVxc 



(12) 



where 
o-ij (r, t) - 



— OO 

2 
3 



V(n(r,t),t-t>) f^pl + ^fn 
V • u{r, t')S i0 



+C(n(r, t),t- t')V ■ u(r, t')8 k 



dt' 



(13) 



77(77., t — t') = f fj(n, u>) exp(—iuj(t — t'))du)/2ir, and simi- 
larly for C. Here n(r, t) and u{f, t) are the time-dependent 
values of the density and velocity field. 

Note that our formula for a xc is still linear in u (f , t) . 
This happens because, due to the constraint (ii) of gen- 
eralized galilean invariance, the velocity must enter the 
stress tensor through its spatial derivatives, which are as- 
sumed to be small, even if the velocity itself is not small. 
Terms of higher order in the velocity would necessarily be 
of higher order in the gradient expansion. For the same 
reason, one can ignore the velocity-dependent terms in 
the xc force (eq. (p!l|)), and the difference between r and 
the "retarded position" R of the fluid clement. Similarly, 
the apparent ambiguity of whether the density entering 
the viscosity coefficients in eq. ([[3]) should be evaluated 
at time t or at some earlier time t', is resolved by noting 
that the difference n(r, t') — n(r, t) = J., V • j (r, r)dr gen- 
erates a higher order gradient correction, provided that 
the range of times which contribute significantly to the 
integral in eq. (13) is essentially finite. 

The simple form of eq. (|l^) is justified by our basic 
assumption that the gradients of the density and veloc- 
ity fields be small. By using the full expression ([ll|) for 
the electromagnetic force on the left hand side of eq. (|T^), 
and by replacing (r, t') — > (R(t'\r, t), t') on the right hand 
side of eq. (|l3|), the approximation can be systematically 
improved, so as to satisfy the local "zero force" theorem 
to all orders in the gradients. A straightforward gener- 
alization allows one to construct an approximation that 
also satisfies the local "zero torque" theorem. 

Finally, we wish to comment on the condition k,q -c 
uj/vp which defines the limit of validity of our approxi- 
mate treatment of a xc . This restriction is forced on us 
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by the analytic structure of the functions fxcLfT)^' 1 ^) 
which are singular along the line u> = kvF, thus limiting 
the radius of convergence of the small-fc expansion. The 
condition is reasonably well satisfied at the characteristic 
frequencies of collective excitations, but it becomes in- 
creasingly restrictive as we lower the frequency, entering 
the domain of single electron-hole excitations. Because 
the region co/kvp — > is not analytically connected to the 
region k — ► 0, it seems unlikely that any local density ap- 
proximation can provide a physically sound description of 
dynamical exchange and correlation in the electron-hole 
excitation region. 

In conclusion, we have developed a consistent and 
physically transparent formulation of nonlinear time- 
dependent density functional theory beyond the ALDA. 
The main physical manifestation of dynamical exchange 
and correlation effects is the appearance of viscosity 
and dynamical shear coefficients in the electron fluid. 
This formulation allows a first-principle calculation of the 
linewidth of high-frequency collective excitations, which 
are not Landau-damped, and would therefore appear as 
sharp (5-functions within the ALDA. As a first appli- 
cation, we have obtained a compact expression for the 
linewidth L of such collective modes, to first order in the 
viscosity coefficients. The result is 

r _ Re J dfu*(f, uj)V k a xc ^ k (f, uj) ^ 
m J drno(r)\u(r, uj)\ 2 

where u(r) is the velocity field of the collective mode cal- 
culated within the ALDA, and no(r) is the equilibrium 
density. Detailed applications of this formula will be re- 
ported elsewhere. 
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